Vibrations of free and embedded anisotropic elastic spheres: 
Application to low-frequency Raman scattering of silicon nanoparticles in silica 
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Vibrational mode frequencies and damping are calculated for an elastic sphere embedded in an 
infinite, homogeneous, isotropic elastic medium. Anisotropic elasticity of the sphere significantly 
shifts the frequencies in comparison to simplified calculations that assume isotropy. New low fre- 
quency Raman light scattering data are presented for silicon spheres grown in a Si02 glass matrix. 
Principal features of the Raman spectrum are not correctly described by a simple model of the 
nanoparticle as a free, isotropic sphere, but require both matrix effects and the anisotropy of the 
silicon to be taken into account. Libration, not vibration, is the dominant mechanism. 
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I. INTRODUCTION 

The classical continuum mechanical problem of the vi- 
brational frequencies of a homogeneous, isotropic, free 
sphere was solved exactly long ago in terms of roots 
of spherical Bessel functionsi. It finds application in 
seismology^, foam spheres^ and even kidney stones^. 
Progress has also been made on generalizations of this 
problem. Closed form solutions are obtainable for spe- 
cial cases of inhomogeneity with certain special kinds 
of anisotropy^, and isotropy with radial inhomogeneity^. 
Generalized anisotropy can be handled for a sphere and 
other shapes with a basis function expansion''. 

Application of the Lamb solution to metal or semicon- 
ductor spheres a few nanometers in size has been suc- 
cessful in explaining vibrational frequencies observable 
using low frequency Raman light scattering^ and other 
methods^. These "nanoparticles" are not truly free, but 
are embedded inside a glass matrix and frequently are 
nearly spherical. The Lamb solution is useful because 
of the sharp change of acoustic impedance z = pc at 
the nanoparticle-glass boundary as long as the nanopar- 
ticle is much heavier or much harder than the glass. The 
applicability of this is challenged in the case of silicon 
nanoparticles grown in Si02 glass, where the densities 
are almost the same and the speed of sound in silicon is 
only 1.5 times faster than in glass. Also, silicon is far 
from isotropic with Zener anisotropy factor 1.6. 

The nanoparticle's vibrational frequencies become 
complex valued due to damping of the modes as energy 
is mechanically radiated away from the nanoparticle into 
the surrounding glass matrix^°. The effect of coupling 
the nanoparticle to the matrix is to shift the frequencies 
of the modes and dampen them.^^. Some qualitatively 
new vibrational modes appear even in the limit of weak 
couplingiSi. The solution to the elastic mechanical prob- 
lem of an elastically anisotropic sphere embedded in an 
infinite elastic continuum has not been presented before. 
Here we employ both a closed form solution for the case of 



isotropic elasticity as well as molecular dynamics (MD) 
simulations for the case of a cubic crystal with elastic 
constants Cn, C12 and C44. The agreement of these two 
methods for the exactly solved free sphere case allow us 
to be confident in the correctness of our solution. Two 
different computer programs were written for the embed- 
ded isotropic sphere calculation. These cross checks are 
essential to avoid errors that have occurred elsewhere in 
the literature of this subject. 



II. THEORY 

The displacement field u{f, t) of an elastic medium of 
density p(r) is governed by Navier's equation: 



CijkljUkd + CijklUkdj = pUi 



(1) 



where Cijki{r) is the 4th rank elastic tensor field with 21 
independent components. In a homogeneous, isotropic 
medium, Eq. Stakes the form of Eq.|21where A and n are 
Lame's constants. 



(A + 2p) V(V ■ u) - fj,y X X uj ^ pu (2) 

To elucidate the key features and physical basis of the 
motion, the vibrations of the nanoparticle are calculated 
here using three methods: (A) analytic solution of a free 
isotropic sphere reprising referencei; (B) analytic solu- 
tion of an isotropic sphere embedded in an infinite ma- 
trix; (C) MD computer simulation of a free sphere includ- 
ing anisotropy. Method (B) has previously been used 
only for some spheroidal modesSii^ or limiting casesiS. 
Method (C) is reported for the first time. 

The following comments apply to methods (A) and 
(B) which are able to exploit spherical symmetry. The 
time dependence of the displacement is taken as u{f, t) = 
■u(r)e*"*, where uj can be complex. The angular fre- 
quency ui (in rad/s) is related to wavenumber i/ (in cm~^) 
through ly ~ cj/(2007rc), where c is the speed of light (in 
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m/s). In what follows, all of the results are reported in 
terms oi Re{v) (energy) and Im{i>) (damping). 

To solve Eq.|21 a scalar and a vector potential are intro- 
duced which correspond to dilatational and equivolumi- 
nal motions respectively. Eq. El separates into scalar and 
vector parts. Detailed expressions for the displacement u 
and the stress tensor can be found in^'*. In a previous 
work£, we expressed u and a • r in a more convenient 
basis. Next, boundary conditions corresponding to the 
system have to be applied and result in secular equations 
whose complex valued roots gives the frequencies of the 
eigenmodes of vibration. Within the nanoparticle, the 
scalar and vector potentials involve the spherical Bessel 
function of the first kind ji . 

For method (A) there is no external force on the sphere 
so that *a-f = at every point on the surface. For method 
(B), the continuity of u and fr-f is imposed at the surface 
of the nanoparticle where it contacts the matrix. 

If the glass matrix is a free sphere of radius i?„i concen- 
tric around the nanoparticle of radius Rp, normal mode 
frequencies can be found^. But what happens is that 
every eigenfrequency vanishes as as Km ^ oo. In 
other words, for a macroscopic glass matrix there is a 
continuum of states but no discrete modes. 

The approach of method (B) is to look for damped 
modes of vibration of the nanoparticle where u) is now 
a complex number such that Re{Lj) is the frequency of 
the vibration and Im{uj) is the damping. In the glass 
matrix, the spherical Hankel function of the second kind, 
hi {kr), is used because it corresponds to an outward 
travelling wave. The wavevector k is also complex. 

Equations for the embedded sphere case of method (B) 
have already been written by Dubrovskiy et al.^". In or- 
der to find the roots of the secular equation, Dubrovskiy 
et al. made approximations which depend on the rela- 
tive values of the parameters in the sphere and around 
it. We made no approximations and found the complex 
roots numericallySS. Frequencies and damping in^« were 
reproduced with this method. 

Method (C) is a molecular dynamics calculationi^. 
While this approach can readily be extended to the more 
general case of an object of arbitrary shape with inhomo- 
geneous anisotropic elasticity, we choose for brevity only 
to present the case of a spherical object with cubic crys- 
talline elasticity. A simple cubic a x a x a lattice of N 
identical point particles of mass pa^ interacting through 
" springs" has particle positions fj integrated in time and 
Fourier transformed in order to obtain mode frequencies. 
Each particle has six first neighbors and twelve second 
neighbors. Every pair (i,j) of first neighbors are cou- 
pled by potential energy ^fcspi(||ri ~ o)^: while sec- 
ond neighbor pairs are coupled by ^ksp2{\\'ri~ rj\ \ — 
•\/2a)^. Cubic octets are coupled by (fc8pt/4a^)(£' — 6a^)^ 

where D = X]i=i \\Rcm - and Rem = |Z]i=i^i- 
The bulk modulus is B ^ (Cn -I- 2Ci2)/3. A smah 
isotropic dilatation S = dV/V stretches all neighbor dis- 
tances proportionally and leads to stored energy ^VBS^, 



so that we obtain a(Cii -I- 2C12) = kspi + 4fcsp2 -I- 24^8^*. 
Considering that a transverse plane wave moves along 
the X-axis at speed \fC^Alp-: we find 06*44 = fcsp2- Fi- 
nally, considering that a longitudinal plane wave moves 
along the x-axis at speed -^/Cn/p, wc find aC\\ — 
kgpi + 2ksp2 + Sfcgpt- Inversion allows the three force 
constants to be obtained for a general cubically elastic 
material: fc^pi = a(Cii — (C12 -I- C44)), ksp2 = aC^ and 
fcspt = |a(Ci2 — C44). 

The frequency dependence of a given mode is f{N). 
A sequence of simulations with N varying from a few 
hundred to 210000 is used to extrapolate /(oo), which 
is the continuum limit. Accurate results are obtainable 
with just a few thousand. Actual nanoparticles are not 
continuum objects, and so the continuum limit is an ide- 
alization. Furthermore, our calculation uses the bulk 
elastic coefficients of silicon. Surface effects (arguably, 
the entire bulk of the nanoparticle is close to its surface) 
and dispersion (effects of the silicon crystal structure on 
phonon speed) are not taken into account in this calcula- 
tion. However, the wavelengths of the lowest modes are 
so large compared to interatomic distances that disper- 
sion will be negligible. 

The correctness and degree of convergence of method 
(C) was checked by comparing its results to those of 
method (A) for isotropic elastic materials of various Pois- 
son ratios. The results of the two methods agree closely. 

Since the first report of low-frequency Raman scatter- 
ing from nanoparticles^, such measurements have been 
made on many materials^. Results are most often 
successfully explained with the free sphere model de- 
scribed above. However, calculated frequencies are much 
too high in silicon nanoparticles embedded in a silica 
matrijiiSiii. In this work, we want to see if the situation 
can be improved by taking into account the presence of 
the matrix and the elastic anisotropy of silicon. 



III. EXPERIMENT 

New experimental results using the samples studied 
iiikL have been obtained. The experimental setup is 
a Jobin-Yvon T64000. Incident light comes from the 
514.5 nm line of an Argon ion laser and is focused onto a 
1 fim^ area on the sample. Power was kept below 1 mW 
to avoid sample heating. The CCD detector allows for 
long accumulation times; the quality of the spectra is 
much improved compared to previous measurements us- 
ing a scanning monochromator coupled with a photo- 
multiplier. Spectra obtained for the sample heat treated 
at 500 °C are presented in Fig. ^ Measurements were 
made with the sample at room temperature. Average 
nanoparticle diameter measured by transmission electron 
microscopy is 6.8 nm. The two narrow lines centered 
around 23 and 29 cm~^ are plasma lines from the laser. 
The parallel and crossed spectra correspond to scattered 
light intensity with polarization parallel (/y) and perpen- 
dicular {I±) to the excitation beam respectively. 
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FIG. 1: Parallel (solid line) and crossed (broken line) config- 
uration Raman spectra of 6.8 nm diameter sificon nanoparti- 
cles embedded in silica. A decomposition is made into peaks 
with depolarization factor 0.75 (plus signs) and 0.33 (bullets). 
Laser plasma lines are indicated by arrows. 

The positions of the peak maxima of these two spectra 
differ slightly indicating scattering by several vibrational 
modes, one being strongly depolarized at low frequency 
and another less depolarized contribution at higher fre- 
quency. From the ratio of the two spectra at 10 cm~^ 
the depolarization factor of the first mode is close to pdep 
= = 0.75, as must be the case for a non-resonant 

non-totally symmetric mode. If we assume our spectra 
components to have only two different depolarization fac- 
tors, then the shape of the scattering with depolarization 
factor different from 0.75 is given by /|| * 0.75 — I±. As 
a result, two vibrational modes are observed: the first 
one with a maximum at 10 cm~^ (determined from the 
crossed configuration) and the second one with a broad 
maximum around 20 cm~^ (determined from the calcu- 
lated spectra). 

IV. DISCUSSION 

First, we want to see how the free sphere model (A) 
compares with our experimental results. Tabled presents 
the calculated vibrational energies for a free isotropic 
sphere of diameter 6.8 nm with the density of silicon and 
different sound speeds corresponding to important crys- 
tal directions in silicon. The first and second [110] lines 
correspond to transverse phonons polarized along [TlO] 
and [001] respectively. Uniformly averaging over all di- 
rections, the longitudinal and transverse speeds of sound 
in silicon are (9017 ± 233) m/s and (5372 ± 385) m/s, 
respectively, where one standard deviation of the varia- 
tion is indicated. The average rms speed is nearly the 
same. For comparison, energies were also calculated us- 
ing these averaged sound speeds. Elastic parameters were 
taken fromi^. As reported by Fujii et al."'^ ^, these frequen- 
cies can't explain the lowest band position. Nevertheless, 



TABLE I; Energy [v in cm^^) of the fundamental spheroidal 
(SPH) and torsional (TOR) vibration mode of a 6.8 nm diam- 
eter free silicon nanoparticle for different sound speeds using 
methods (A) and (C). Im{v)=0. 



Calculation 


Speed of 




SPH 




TOR 


method 


sound 


/=0 


1=1 


1=2 


1=1 1=2 


(A) 


[100] 


28.3 


26.4 


23.6 


52.5 22.8 


(A) 


[111] 


38.7 


28.0 


21.1 


45.8 19.9 


(A) 


[110] 


38.8 


26.0 


19.3 


42.0 18.3 


(A) 


[110] 


33.8 


29.1 


23.7 


52.5 22.8 


(A) 


average 


35.4 


28.2 


22.0 


48.3 21.0 


(C) 


anisotropic 


34.9 


27.0 


19.4 


50.4 19.4 



we will use these results to choose the isotropic parame- 
ters which best describe the calculated anisotropic ener- 
gies. From Tabled the closest agreement with the exact 
anisotropic calculation (C) is obtained when method (A) 
uses the averaged sound speeds. 

We applied method (B), with averaged sound speeds, 
to silicon nanoparticles embedded in silica glass. Elastic 
parameters for silica were taken fromi^ andiS. Results 
are presented in Table ITTI Results using method (A) (free 
sphere) are also shown, for which case there is no damp- 
ing. Modes can be spheroidal (SPH) or torsional (TOR). 
I is the orbital quantum number and n is the index of the 
solution. Blank spaces in the free case indicate that there 
is no equivalent to the corresponding embedded mode. 

The silica matrix has three effects: it shifts ener- 
gies (i?e(i/)), it creates damping {Im{v)) and it also 
introduces new modes. One of these new modes is 
(TOR,/=l,n=0). In the absence of the matrix this would 
correspond to rigid rotation. The matrix imposes a 
restoring force changing the rotation into a libration 
(rigid rotational oscillation). Other new modes have 
more complicated origins. Considering modes existing 
in the free case, the more their energies are shifted, the 
greater the damping. This is consistent with a stronger 
contribution of the matrix. The damping values given 
in Table [H] correspond to the Raman homogeneous half 
width at half maximum. The observed Raman linewidth 
should be greater mainly because of the nanoparticle size 
distribution inside our samples. 

For comparison purposes, calculated energies and 
damping of vibrations of a 6.8 nm diameter spherical cav- 
ity in silica as well as around an infinitely hard, heavy 
sphere are presented in Table lllll In particular, the 
lowest spheroidal modes in Table ITU correspond approxi- 
mately to modes around a rigid, heavy sphere. 

From Table ^1 two modes are good candidates 
for the lowest Raman peak: (SPH,^=l,7i=l) and 
(TOR,;=l,7i=0). But SPH modes with odd I and TOR 
modes with even I have odd parity and should be Raman 
inactive^. The calculated damping for (TOR,Z=l,n=0) 
is higher than the Raman peak HWHM. Moreover, size 
and shape distributions^- and temperature broadening 
also increase the Raman linewidth. A possible explana- 
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TABLE 11: Energy {Re{u) in cm"^) for free (method (A)) 
and silica embedded (method (B)) 6.8 nm diameter sihcon 
nanoparticles for various modes. In the embedded case damp- 
ing {Im[u) in cm~^) is the half width at half maximum. 



mode 




7 
( 


n 




free (A) 
energy 


embedded (B) 
energy damping 


SPH 








1 
2 




35.4 
85.1 
130.5 


38.3 
85.5 
130.8 


13.9 
10.9 
10.5 






1 



1 
2 




28.2 


6.3 
9.3 
32.7 


9.9 
16.2 
11.4 






2 



1 
2 




22.0 
40.2 


18.8 
22.5 
46.4 


18.2 
10.1 
14.0 


TOR 




1 



1 
2 




48.3 
76.2 


11.8 
48.4 
76.3 


11.0 
6.8 

6.7 






2 



1 
2 




21.0 
59.7 


6.5 
26.1 
59.9 


12.2 
13.7 
7.3 


TABLE III: Energy and damping of a silica matrix surround- 
ing a 6.8 nm diameter spherical cavity and an infinitely rigid 
and heavy sphere. Only even parity modes are shown. 


mode 


I 


n 


cavity 
energy damping 


rigid heavy sphere 
energy damping 


SPH 










9.1 


7.4 








2 



1 
2 


6.1 
7.8 
23.3 


4.6 
20.3 
11.1 


15.2 


14.5 


TOR 


1 







5.1 


8.8 







tion is that in our model, the nanoparticle-matrix contact 
is perfect and therefore more energy is radiated into the 
matrix than in our sample. 

The second mode is in agreement with the calculated 
(SPH,Z=2,n=0,l) modes. FoUowingSi, we assume the de- 



polarization factor of this peak to be 1/3. Therefore, we 
are able to fully decompose each spectra. The decompo- 
sition for the parallel configuration is shown in Fig. ^ 



V. CONCLUSION 

The special demands of silicon nanoparticles embedded 
in silica require consideration of mode broadening by the 
matrix, mode shifting by the matrix, new modes created 
by the matrix and the anisotropy of the silicon. A strong 
candidate as the dominant feature of the Raman spectra 
is the libration mode (TOR,Z=l,n=0) which would be 
Raman inactive without either anisotropy or nonspheri- 
cal shaped. 

The approach we are using can easily be extended to 
handle other features of realistic nanoparticles, such as 
non-spherical shape, surface relaxation and details of the 
nanoparticle-matrix interface. The extension to inhomo- 
geneous onion-like particles where the parameters have 
radial dependence can also be made. 

Exact incorporation of anisotropy is not difficult, but 
relatively accurate estimates of frequency are possible us- 
ing an isotropic model of the elasticity as long as direction 
averaged speeds of sound are used. This works best for 
modes with I <\, since they are at most three- fold degen- 
erate with isotropy and are not split by cubic crystalline 
anisotropy. 

For an isotropic material, the 1—2 spheroidal and tor- 
sional modes are 5-fold degenerate. This degeneracy is 
broken by the anisotropy, into what MD simulations sug- 
gest is a doublet. 

Overall, this study has provided a clearer basis for 
interpreting Raman experiments where the acoustic 
impedance of the nanoparticle is close to that of the 
matrix. Crystal elasticity is anisotropic in any case. 
However, without theoretical calculations of the abso- 
lute cross sections of Raman scattering it is possible only 
to make tentative identification of the vibrational modes 
that appear in the low frequency Raman spectrum. 
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